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ABSTRACT 

A  gas -cooled,   cylindrical,   nuclear  reactor  is  used  as  the 
basis  in  the  theoretical  derivation  of  the  coolant's  enthalpy  increase 
and  pressure  loss  across  the  reactor.     Turbulent  flow  is  assumed  in 
the  coolant  passages  of  the  reactor,   and  the  Reynolds  analogy  is 
used  as  the  basis  for  correlating  the  heat  transfer  coefficient  and 
the  friction  coefficient. 

The  general  equations  are  derived  and  two  examples,   a 
nuclear -hydrogen  rocket  and  a  ramjet,   are  given  to  demonstrate 
applications  of  the  general  results. 


SUMMARY 

The  development  of  nuclear  reactors  has  led  to  a  new  class  of  heat 
transfer  problems.     Heat  is  produced  within  a  solid  fuel  at  a  rate  that 
is  almost  independent  of  the  coolant  flow  rate,   at  least  for  gas -cooled 
reactors.     The  power  density  varies  with  position  in  general,   although 
it  can  be  made  nearly  uniform  with  the  introduction  of  appropriate  re- 
flectors and  with  variation  of  fuel  concentration.     However,   these 
measures  are  not  desirable,  particularly  if  there  are  weight  restric- 
tions for  the  reflector  and  cost  restrictions  on  fuel  elements. 

Aside  from  considerations  of  weight  and  cost,   the  most  important 
characteristic  of  a  reactor  is  the  ratio  of  pumping  work  to  the  increase 
of  thermodynamic  availability  of  the  coolant  passing  through  the  re- 
actor.    For  any  single  coolant  passage  through  the  reactor,  this  ratio 
can  be  determined  readily  from  the  Reynolds    analogy,  for  any  speci- 
fied power  distribution  along  the  length  of  the  channel,   and  a  given 
coolant  flow  rate  through  the  channel  or  a  given  pressure  drop  over  the 
reactor.     The  limitation  on  the  solid  fuel  temperature  will  fix  the 
minimum  pressure  drop  or  flow  rate.     If  the  power  density  is  the  same 
for  all  channels,   the  results  are  simply  summed  over  the  reactor 
cross-section. 

When  the  power  density  varies  from  channel  to  channel,  there  is 
no  simple  way  of  calculating  pressure  drop  or  temperature  rise 
through  the  reactor.     The  pressure  drop  for  all  channels  must  be  the 
same,   but  the  flow  rate  and  temperature  rise  will  vary  from  channel  to 
channel.    Apparently  the  calculations  of  temperature  and  flow  rate  dis- 
tributions over  the  reactor  have  been  done  by  numerical  computation 
in  the  past.     An  iterative  procedure  is  developed  in  this  thesis  for 


determination  of  the  overall  reactor  performance  in  analytical  form. 
The  sole  restriction  on  the  procedure  ie  that  the  Mach  number  of  the 
flow  entering  any  coolant  channel  be  appreciably  less  than  unity.     This 
restriction  is  not  a  serious  one  because  the  pressure  drop  resulting 
from  higher  entering  Mach  numbers  would  be  generally  unacceptable. 

The  principal  objective  of  this  thesis  is  the  development  of  the 
method  of  calculating  reactor  performance  for  any  specified  power  den- 
sity distribution.     In  final  form  the  method  gives  the  pressure  drop, 
the  temperature  distribution  and  mass  flow  rate  distribution  in  the 
coolant  leaving  the  reactor,   as  well  as  maximum  wall  temperatures, 
where  porosity,   coolant  channel  geometry,   and  power  density  distribu- 
tion are  prescribed.     The  temperature  of  the  coolant  leaving  the  re- 
actor can  be  made  uniform  over  the  reactor  cross- section  by  introducing 
appropriate  restrictions  in  the  individual  coolant  channels.     This  pro- 
cedure increases  the  average  temperature  rise  through  the  reactor  for 
a  given  maximum  wall  temperature,   but  does  so  at  the  cost  of  in- 
creased pressure  drop.    Whether  restriction  of  the  coolant  flow  channels 
in  this  way  is  advantageous  or  not  will  depend  upon  the  particular  ap- 
plication.    The  method  of  analysis  developed  here  can  be  used  to  deter- 
mine the  optimum  restriction  of  the  flow  channels  for  any  specific 
example. 

To  illustrate  the  application  of  the  method  for  calculating  reactor 
performance,   the  nuclear  rocket  and  nuclear  ramjet  were  chosen  as 
examples.     The  specific  impulse  for  the  rocket  is  found  to  be  higher 
when  the  coolant  passages  are  restricted  to  give  a  uniform  coolant  out- 
let temperature  than  for  no  restriction.     The  thrust  coefficient  of  the 
ramjet  is  not  improved  appreciably  by  restricting  the  flow  rate. 
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SYMBOLS 

A  Coolant  passage  cross -sectional  area 

&-  Radius  of  reactor 

C,  Friction  coefficient 

C,  Heat  transfer  coefficient  (Stanton  number) 

c  Specific  heat  at  constant  pressure 

f(£)  Axial  power  distribution  function 

g(  rj  )  Radial  power  distribution  function 

k  Conductivity 

Coolant  passage  length 

Total  mass  flow  rate  through  reactor 

rh  Mass  flow  rate  per  unit  area 

P  Coolant  passage  wetted  perimeter 

p  Gas  pressure  (static) 

q  Power  density  i.i  reactor;  heat  produced  per  unit  volume 

q  Maximum  power  density  in  reactor 

R  Gas  constant 

T  Gas  temperature  (static,   averaged  over  coolant  passage 

cross-section) 

u  Gas  velocity  (averaged  over  coolant  passage  cross - 

section) 

^  Porosity;  void  volume  /total  volume 


I 

rh 


p/ 

C,  — y-  i   heat  transfer  parameter 


p,  f    F/    ,   friction  parameter 

7  r/a.    ,   dimensionless  radius  from  centerline  of  reactor 

/  Viscosity 

^  Kinematic  viscosity 

£  x//    ,   dimensionless  distance  from  coolant  passage 

entrance 


SYMBOLS 

p  Gas  density  (averaged  over  coolant  passage  cross 

section) 

&~  Prandtl  number 

Subscripts: 

0  Entrance  to  coolant  passage 

1  Exit  from  coolant  passage 
w  Wall  of  coolant  passage 

T  Total 


I.     INTRODUCTION 

The  power  distribution  in  a  nuclear  reactor  is  proportional  to 
the  neutron  density  which,   in  a  bare  cylindrical  reactor,   is  highest  in 
the  center  and  decreases  both  axially  and  radially  away  from  the  cen- 
ter.    When  the  reactor  is  used  as  a  source  of  heat,  the  variation  in 
power  density  distribution  results  in  an  uneven  heat  supply  to  the 
coolant. 

The  neutron  density  can  be  made  more  uniform  by  enclosing  the 
reactor  with  a  reflector  which  causes  some  of  the  leaking  neutrons  to 
diffuse  back  into  the  core,   or  by  varying  the  fuel  loading  so  that  the 
fuel  concentration  is  higher  towards  the  edges  than  in  the  center. 
However,  both  of  these  methods  have  disadvantages.     The  reflector 
material  is  usually  quite  bulky  and  relatively  heavy.     As  a  result,   in 
several  applications  of  nuclear  reactors  as  a  source  of  power,   spe- 
cifically,  for  airborne  propulsive  applications,  it  is  necessary  to 
restrict  the  amount  of  reflector  material  that  might  normally  be  used. 
Fuel  concentration  variation  is  undesirable  since  reactors  are  usually 
constructed  by  assembling  many  small  lamina,   or  rods,  to  obtain  the 
final  configuration.     For  manufacturing  ease  and  interchange  ability, 
it  is  more  desirable  to  have  'uniformity  of  the  fuel  elements. 

The  uneven  power  density  distribution  of  the  nuclear  reactor 
results  in  the  coolant  having  temperature  gradients  in  tlie  exhaust  if 
the  mass  flow  in  each  of  the  coolant  channels  is  the  same.     The  cool- 
ant in  the  channel  through  the  center  of  the  core  will  be  heated  the 
most,  whereas  the  flow  in  the  outer  channels  will  be  much  cooler  due 
to  the  lower  heat  generation  at  the  reactor  extremities.     Flow  of  this 
type  is  usually  undesirable,    since  the  surface  temperature  of  the  chan- 
nel walls  in  the  outer  channels  is  much  lower  than  the  maximum 
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allowable,   and  thus  the  overall  efficiency  of  the  power  plant  is  reduced. 
In  addition,   the  temperature  gradients  of  the  exhaust  of  this  type  flow 
lead  to  mixing  losses  in  the  nozzle. 

Theoretically,   these  exhaust  gas  temperature  gradients  can  be 
eliminated  by  restricting  the  mass  flow  in  the  outer  channels  so  that 
the  mass  flow  in  every  channel  is  proportional  to  the  heat  release  in 
that  channel.     The  result  is  that  the  outer  channels,   where  the  heat 
release  is  less,   are  cooled  to  a  lesser  degree  so  that  each  channel 
operates  at  its  maximum  allowable  wall  surface  temperature,   and  the 
temperature  of  the  exhaust  gases  from  the  reactor  is  uniform.     This 
mass  flow  restriction  can  be  accomplished  by  varying  the  flow  channel 
cross- sectional  areas  continuously  across  the  reactor  core  so  that  for 
the  same  overall  pressure  drop  across  the  channel  the  mass  flow  is 
smaller.     The  same  effect  can  be  obtained  by  installing  orifices  in  the 
channels  to  induce  a  pressure  drop  so  that  again  the  mass  flow  is  pro- 
portional to  the  heat  release  of  the  particular  channel.     Of  these  two 
methods,   the  former  is  preferable  since  it  utilizes  all  the  available 
pressure  drop  to  improve  the  heat  transfer  coefficient.     However,   it 
would  not  usually  be  practical  to  assemble  a  reactor  in  this  manner, 
and  so  the  orificing  technique  is  considered  to  be  more  practical. 
Elimination  of  the  exit  temperature  gradients  and  increasing  the  wall 
surface  temperatures  of  all  channels  to  the  maximum  allowable  by 
either  of  these  methods  increases  the  overall  pressure  drop  across 
the  reactor  and  reduces  the  overall  mass  flow,   in  addition  to  increas- 
ing the  average  temperature  rise  per  unit  mass  flow.     Because  of  the 
increased  pressure  losses  and  the  reduction  in  mass  flow,   the  over- 
all performance  of  the  reactor  must  be  evaluated  to  determine  the 
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performance  gain,   if  any,   by  orificing  tc  ascertain  if  the  increased 
design  and  assembly  complications  are  warranted. 

W.   B.    Hall  (1)    has  worked  out  an  example  of  the  restricted  flow 
reactor  configuration.     The  basic  reactor  assumptions  made  by  Hall 
and  those  of  this  the3i3  are  similar.     However,   the  analysis  here  is 
more  general,   including  allowance  for  higher  pressure  drop  and 
higher  velocities  in  the  coolant  channels,   arid  a  procedure  for  finding 
overall  performance  when  the  coolant  outlet  temperature  varies 
radially.     These  improvements  in  the  analysis  are  needed  for  appli- 
cation to  reactors  for  propulsion. 

In  order  to  keep  the  general  discussion  of  heat  transfer  and 
pressure  drop  in  gas -cooled  reactors  for  propulsive  applications  with- 
in reasonable  limits,   simplifying  assumptions  are  required.     For 
this    report,  the  reactor  was  taken  a3  cylindrical,  with  radial  sym- 
metry in  fuel  loading  and  coolant  passage  distribution.     The  coolant 
passages  were  assumed  to  be  of  uniform  size,  and  the  cross -sectional 
area  was  assumed  constant  over  the  length  of  the  reactor.     The  power 
density  of  the  reactor  was  assumed  to  vary  axially  and  radially  ac- 
cording to  simple  algebraic  laws.     Since  the  analysis  was  made  pri- 
marily to  investigate  the  result  of  using  a  reactor  as  a  type  of  heat 
exchanger,   or  source  of  heat,   reactor  physics  as  such  was  not  con- 
sidered,  and  no  account  was  taken  of  the  power  distribution  shift 
during  operation,   etc.     The  porosity  of  the  reactor  was  not  assumed 
constant  in  the  derivation  of  the  general  equations,   but  was  assumed 
constant  in  specific  examples  worked  out  to  demonstrate  the  use  of 
the  results. 
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The  evaluation  of  the  overall  performance  of  a  propulsion 
system  involves  determining  the  overall  enthalpy  rise  of  the  working 
fluid  and  the  losses  of  the  system.     For  a  reactor  of  given  design, 
the  controlling  factor  in  operating  power  is  normally  some  maximum 
allowable  temperature  in  the  system.     ThiG  limit  may  be  set  by  the 
phase  stability  of  the  fuel  elements,   or  the  moderator,   by  the  allow- 
able thermal  stresses  in  the  fuel  or  some  other  part  of  the  system,   by 
the  influence  of  temperature  on  corrosion,    or  some  other  particular 
thermal  effect.     For  a  given  design  and  this  limiting  temperature, 
the  temperature  distribution  in  the  reactor  core  can  be  determined, 
and  from  thi3,   the  maximum  allowable  coolant  wall  temperature 
subject  to  the  limiting  temperature  of  the  particular  configuration. 
Since  this  report  was  intended  to  be  quite  general,   the  controlling 
factor  for  reactor  operation  was  selected  as  the  maximum  allowable 
coolant  wall  temperature  which  can  be  determined  when  given  a  par- 
ticular configuration  with  type  of  fuel,   concentration,   etc.     Thus  the 
problem  to  be  solved  involves  calculating  the  overall  enthalpy  rise 
and  preesure  losses  for  a  given  maximum  wall  temperature  and  given 
initial  conditions. 

The  high  power  requirements  of  today's  rockets  and  ramjets 
and  the  expanding  technology  in  the  nuclear  reactor  field  make  the  pos- 
sibility of  using  the  gas -cooled  reactor  as  a  source  of  propulsive 
power  for  these  vehicles  a  distinct  reality.    Applications  for  rocket 
engines  hold  promise,   since  the  effective  exhaust  velocity  is  propor- 
tional to  the  inverse  of  the  square  root  of  the  molecular  weight,   which 
suggests  the  possibility  of  using  low  molecular  weight  gases  such  as 
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hydrogen  for  the  working  fluid.     For  ramjet  applications,  the  nuclear 
reactor  offers  the  possibility  of  relatively  unlimited  range  or  endur- 
ance at  low  as  well  as  high  altitudes. 

In  this  thesis,  the  general  equations  are  derived;  and  two  exam- 
ples,  a  rocket  using  hydrogen  ax  the  working  fluid  and  a  ramjet,   are 
given.     The  results  of  these  examples  show  a  distinct  advantage  in  ori- 
ficing  to  increase  the  specific  impulse  of  a  rocket,   but  the  advantages 
for  a  ramjet  are  so  slight  that  it  is  doubtful  if  the  additional  complexity 
in  reactor  construction  arising  from  orificing  is  warranted. 
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H.     ANA  LYSIS  Aim  GENERAL  SOLUTION 

2.  1    Refolds  Analogy 

In  most  nuclear  reactors  the  flow  of  the  coolant  through    the 
cooling  passages  is  turbulent.     For  this  reason  the  motion  of  the  fluid 
cannot  be  defined  exactly,    so  exact  analytical  solution  of  flow  and  heat 
transfer  problems  is  impossible.     However,   using  the  equations  of 
motion,   where  the  laminar  transport  properties  are  replaced  by  ef- 
fective values  which  account  for  the  effects  of  turbulence  and  which 
must  be  determined  ej-perimentally,  we  are  able  to  solve  these  prob- 
lems. 

The  Reynolds  analogy,   first  postulated  by  Osborne  Reynolds  in 
1890  and  verified  by  experiment,    states  that  momentum  and  energy 
are  transferred  in  the  same  way  in  turbulent  shear  flow.     This  anal- 
ogy can  be  expressed  quantitatively  by  the  equation: 

(2    1) 

where    q      is  the  heat  flow  rate  per  unit  area  normal  to  the  surface 
and     T0      is  the  shear  stress  at  the  wall.     For  fully  developed  turbu- 
lent flow,   the  velocity  profile  normal  to  the  surface  is  nearly  flat; 
hence,   integrating  between  the  wall  where    u    is  zero  and  a  station 
where    u    and   T    can  be  considered  to  be  mean  values 


oCpiT-T^ju.  P 


u 


T  (3.  2) 


or 


CK     =    h   Cf  (2.  3, 
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where    C,     and    Cf   are  the  heat  transfer  coefficient  (Stanton  number) 
and  friction  coefficient  respectively. 

If  we  further  write  the  equations  of  momentum  and  energy, 
averaged  with  respect  to  time,  for  mean  flow  parallel  to  a  surface  we 
find: 

V0    -    -  f>   "V'      +^  du/dy  (2.4) 


fr*     -fCp'V   +    kaT/^  (2.5) 


where   U.  is  the  viscosity  and   k    is  the  conductivity,   and  upon  re- 
arrangement to  better    show  the  similarity: 


-r  (2.6) 


r  p  K   dt/a   '  "^'^j 


where     7)  =  /^    is  the  kinematic  viscosity  and     (T  -/*&  *s  &**  Prandtl 
number.     Thus,   it  can  be  seen  that  if  the  Reynolds  analogy  holds  for 
the  turbulent  properties,   and  if  cr  =  1,   Equation  2.  3  holds  through- 
out turbulent  and  laminar  regions. 

For  most  gases      (T  •— ->    0.  7  -  1.0,   so  the  Reynolds  analogy  can 
be  assumed  to  hold  reasonably  well.     The  analysis  and  examples  be- 
low are  restricted  to  gas -cooled  reactors  with  turbulent  flow  in  the 
coolant  channels. 
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Z.  2    Temperature  Distribution 

Duo  to  non- uniformity  of  the  heat  generation  in  the  core  of  the 
reactor,   the  average  temperature  rise  of  the  coolant  through  the 
reactor  core  must  be  determined  to  find  the  overall  enthalpy  rise 
through  the  reactor.     The  temperature  rise  in  the  individual  channels 
must  be  found  first.     Integration  of  each  coolant  channel  temperature 
rise  over  the  reactor  cross-section  gives  the  total  enthalpy  rise. 

The  heat  balance  for  a  length  of  coolant  passage,     dx  ■  /d£,,   at 

a  given  radial  station,    Y)    ,   is  given  by: 

( 

rncpAdT    =     Jj.Aif(S)  J^-ld5  <2'7> 

where    m    is  the  mass  flow  rate  per  unit  area,     c      is  the  specific 
heat  at  constant  pressure,     A    is  the  cross -sectional  area  of  the 
coolant  passage,     T    is  the  static  temperature  of  the  coolant,     q      is 
the  power  density  at  the  point  of  maximum  heat  generation  in  the 
reactor  core  (usually  at  the  center  of  the  core),   £  =  j   is  the  dimen- 
sionless  coolant  passage  length  measured  from  the  channel  entrance, 
yj  s    —is  the  dimensionless  radius  of  the  core  measured  from  the  cen- 
ter,    f(£)  is  the  power  distribution  function  in  the  axial  direction,    g(rt  ) 
is  the  function  in  the  radial  direction,   and  oi    is  the  porosity,    or  void 
fraction,   defined  as  the  ratio  of  the  total  coolant  passage  cross -section- 
al area  to  the  total  cross -sectional  area.     In  general,    K     could  be  a 
function  of    >?   . 

Rewrite  Equation  2.  7  as: 

*I.Alf($)  (2.7a) 
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where      A  =        n°  -zr    -L— ^     is  a  dimensionless  parameter.     Since    m 

is  constant  for  a  given  coolant  passage.    Equation  2.  7a  can  be  integrated 


to: 

T 


=  Ta|i    +  X(>,)f  f($)d^  (2.8) 


The  heat  balance  between  the  wall  and  the  coolant  is  given  by 

the  relation: 

ny  cp  AiT    -   CK  f>  a.  cp  (T„  -Tj  PJ  d£  (2.9) 

where    C,     is  the  Stanton  number  and  depends  chiefly  upon  the  Rey- 
nolds number  for  gases  with  Prandtl  number  close  to  unity;    p    is  the 
coolant  density,    u  is  the  velocity  of  the  coolant,    T    is  the  channel 
sur face  temperature,   and    P    is  the  wetted  perimeter. 

An  empirical  relation  used  in  flow  through  pipes  of  circular 

cross -section  is: 

/«  o.oz.3 

CK    "      ~^~l      oc  (2-  10) 

where  R      is  the  Reynolds  number.     The     Reynolds  number  is  de- 
fined as: 

P     =    4-^  (2.11) 

where  4A/P  a  D,    ,  the  hydraulic  diameter.     Since    R      for  a  constant 
area  channel  varies  only  with  u.  ,  which  in  turn  depends  only  on  the 
temperature  and  is  approximately  proportional  to  the  three- fourths 
power  of  the  absolute  temperature,   the  variation  of    C,     (and  of  C,  , 
the  friction  coefficient)  along  the  passage  will  be  quite  small,  partic- 
ularly for  turbulent  flow. 
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Since    m  »  pu  ,    Equation  2.  9  can  be  written  as: 

4£    ■    iflfT^-T)  (2.  12) 

P 
where    psC.  -t-    ,   is  constant,   and  is  a  dimensionless  parameter 

proportional  to  the  heat  transfer  coefficient,   C,  .     Hence,   the  wall 

temperature  is  given  by: 


or,   from  Equation  2.  8: 

Tw    -    T0    +AT.JJL£(£)   *    [-f(*)d^  (2.H) 

* 

dTw 
The  maximum  wall  temperature  is  given  by  the  condition,  -j* —  *  o, 

i.  e.  , 

dk  L*^  *  ^  S)  =°  (2#  14a) 

From  this  equation,  the  value  of  £  for  the  maximum  wall  temperature 
can  be  found,   and  substituting  into  Equation  2.  14    Tw      ,   the  maximum 
wall  temperaturejcan  be  obtained. 


w 
m 


2'  ^    Pressure  Distribution 

The  coolant  will  undergo  a  pressure  drop  as  it  flows  through 
the  channels.     This  pressure  drop  arises  from  two  factors;  the  ac- 
celeration of  the  coolant  as  it  passes  through  the  channels,   and  the 
friction  drag  opposing  the  flow  of  the  coolant  through  the  passages. 
The  acceleration  of  the  flow  results  from  decreasing  density  of  the 
coolant,  which  in  turn  results  from  increasing  temperature  and  de- 
creasing pressure.     The  friction  drag  cannot  be  avoided  since  the 
Reynolds  analogy  holds. 
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The  pressure  drop  along  an  individual  channel  can  be  determined 
by  applying  the  equations  of  momentum  and  continuity.  The  momentum 
balance  for  an  element  of  length  of  the  coolant  passage  is  given  by: 

where    p    is  the  static  pressure  of  the  coolant  and    C,   is  an  empiri- 
cal friction  coefficient.     From  the  continuity  and  perfect  gas  relations: 

where    R    is  the  gas  constant  and  subscripts,    o,    refer  to  the  entrance 
station  of  the  coolant  channel.     Substituting  Equation  *"  i&  into  Equa- 
tion *~  15  and  rearranging: 

where    (3,  »   -w—    -y—  ,   a  dimensionlese  parameter  proj>ortional  to  C*, 
the  friction  coefficient. 

The  second  term  in  tlie  coefficient  of    dp/d£    in  Equation    2«  17 
is  a  result  of  the  acceleration  due  to  the  pressure  drop  in  the  channel. 
In  many  applications  this  term  can  be  neglected  and  then  Equation  z*  * ' 
can  be  integrated  exactly:  however,   hi  flow  through  a  reactor  for  pro- 
pulsion applications,   the  pressure  drop  may  be  quite  large,   so  this 
term  must  be  retained.     The  term  containing    dT/d|    is  a  result  of 
acceleration  due  to  the  heat  input  to  the  coolant  as  it  flows  through 
the  passage,   and  will  be  quite  large  for  nuclear  reactors.     Equation 
*■  *■ «    is  an  exact  equation  for  the  pressure  drop,   but  it  cannot  be 
integrated  in  closed  form.     However,   the  equation  can  be  integrated 
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by  an  iteration  process  if  the  Mach  number  is  low. 

2.  3.  1    Iteration  Procedure 

Introduce  a  dimensionless  parameter,  £    ,   defined  as: 

6    =  (mf&J*    -    (p.u.)1  RT0       =      u.^  - 

P.1  Cp.RT.)*        "r^-      ^  l 

The  subscripts,    o,   refer  to  the  channel  entrance,   and  since  the 
product    RT      is  proportional  to  the  square  of  the  local  speed  of 
sound  for  perfect  gases,  £    is  seen  to  be    -y^    »  **"*  thus  propor- 
tional to  the  Mach  number  at  the  channel  entrance.     It  will  be  as- 


sumed that    -Jfi—   =  ( -CM  JJue     «  i      as  well. 

RT  V    pv  )     T 


r        ^  P^T 
Equation  2.  17  becomes,  upon  substitution  of  a  , 

f-f  -  e  J&-1  fcliW£-^(4Ii  (5.  T)  =  o 
Let   p   be  expanded  in  a  power  series  in    e  • 


(2.  18) 


(2.  19) 


Substituting  this  series  into  Equation  2.  18,  we  obtain: 


and,   upon  collecting  powers  of   €.  ,  the  differential  equation  for  p*   ' 
is  found  as: 


&"+M%  +  *r> 


jL  (±  \       J-./  dr  -  a.  r\  =Q  <2-21> 
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and  integrating  Equation  2.  21, 

For  the  second  approximation,  the  equation  for    p*       ie: 


and  integrating  Equation  ?..  23 
rx)  A 


o 

Define  the  static  pressure  loss  as: 


(2.  24) 


e    p.    +  • 


Ap«  K  -  (fb*  feF»   ^  e    P 


(<i 


where     ^>(     is  the  first  approximation  for    p    at  station  1,  the  reac- 


<"») 


tor  exit,   and     pv     is  the  second  approximation  at  this  station.     There- 
fore,  the  pressure  loss  ratio  may  be  expressed  as: 

Aft.    -  *J$L   -.*_£•_ 

P-  P-  IS, 


=    +  <EX  jC*r  ,.  /3+T)«i§  (2.25) 

f  eHV^T){l'^^r'Uf'','r(r)l 4J 


±e.  .  «  f;  +  •»  f, 


(2.26) 


These  functions,    F.    and    F?  ,   may  be  expressed  in  a  more  con- 
venient form  by  integration  and  rearrangement  to  the  following: 
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-?t,(t1-u).3,Az[i;  -  £f^a*] 


(2.  27) 


(2.  28) 


2.  4    First  Order  Approximation 

The  function,     f{|)»  has  been  defined  as  the  axial  power  density 
distribution.     Hence,     F,    may  be  written  as: 

F.  =  h   +  A  L  B°  +  ^  ( s- "  B<  jJ  (2- 29) 


where    B      and    B,    are  defined  as: 
o  1    . 


3.  s    £f(5M?    =   ^"*r- 


(2.30) 


(2.31) 


If  we  further  require  that    j;{£)  be  symmetrical  about  the  center  of 
the  reactor,     £  »  \t  integration  of  3,    shows  that  B.  ■  -*B     .     There- 
fore,    F,  ,   under  these  assumptions  becomes: 

f     =,      ;\     «■  X  B.  (  I  +  /%  )  {2.  32) 


In  order  to  integrate  Equation  2.  22  to  find  the  first  order  ap- 
proximation to  the  static  pressure  loss  across  the  reactor.    Equation 
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2.  7  can  be  rewritten  as : 

%   -   AT.  -f(§)  =  K  C  #£M  -U-    f  I  5)  (2-  33) 

a  5  °^      ye 

/  /  /  _  ,\ 
where    K    is  defined  as:       |K    »     — ^ : — — 

c  |p,  v  i- i; 

and  tlie  parameter,    X    ,  becomes: 

The  exit  temperature  of  the  working  fluid  as  a  function  of  radial 
position  can  be  determined  by  integrating  Equation  2.  33,   since  the 
integration  is  in  the  axial  direction  only: 

Tfy)  -T.  -  KTo  _^     '      B,  <2. 34) 

Using  these  results.   Equation  2.  26  for  the  first  approximation  to  the 
pressure  loss  can  be  expressed  in  either  of  the  following  equations: 

P- 
or,   alternatively,  using  the  relation  of  Equation  2.  34: 

Ae.,*^    +   TiilLJl(i  +/v)«   «-•••■       (2.36) 

r-  t.  y 

Neither  of  the  two  Equations  2.  35  and  2.  36  is  in  a  convenient  form 
for  integration.     However,  without  changing  the  value  of  the  equation, 

we  may  add  and  subtract  the  term,    e.  — '  x     — (. '  *  ftVx.)  ,  from  the 

T     —  T 
right  side  of  the  equation.        — ! —      is  at  present  unknown  and  is 

'  • 
to  be  determined.     Performing  this  operation.  Equation  2.  35  becomes 

"pT  =  *[ff  +  ^*4**ri{  t(^^/j[K6„|^-t^e]    (2.37) 
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For  a  given  reactor  under  given  flow  conditions,       — — — —   will  be  a 
constant.     Further,   when  integrated  over  the  radial  area  of  the  reac- 
tor,    T.    and  T.    will  not  differ  to  a  great  degree,   so  the  second 
term  on  the  right  side  of  Equation  2.  37  will  now  be  very  small  when 
compared  to  •Uher  of  the  other  two  terms. 

For  convenience,   rewrite  Equation  2.  37  as: 


p° 


(2,38) 


or,   rearranging  in  terms  of     y6      : 


c&)\[±&-%(Grj 


-i  -  *?/? 


—     -o 


(2.39) 


where: 


(2.  40) 


Since  the  middle  term  in  Equation  2.  39  has  been  purposely  made 
very  small  in  comparison  with  the  other  two  terms,   this  equation  can 
now  be  solved  by  successive  approximations.     Write  Equation  2.  39  as: 


lm 


.[(Br^K*"-^^"^ 


=  o 


(2.41) 


and  assume  a  solution  of  the  form: 

/r  ■  '•"  r 


a. 


(,  * 1.5"  +y>^+ )"| 


(2.  42) 


After  substituting  this  assumed  solution  into  Equation  2.  41,   col- 
lecting similar  powers  of    &    ,  and  taking  the  limit  as  J"->.  I  ,   we 


obtain: 


■l—y 
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c  D 


P  (2.43) 


+  t 


C  S  CD 

i-    3 


_EAP/po    £/e^/Po 


It  is  more  compact  to  leave  Equation  2.  43  in  terms  of  C,   D,   and  E, 
since  they  are  defined  in  Equation  2.  40.     For  integration  in  the  radial 
direction,   it  can  be  :seen  that  C  contains  the  variable  of  integration, 
g(  V|  ).     For  the  first  approximation,   use  the  first  term  (which  does 
not  contain  C  or  D);  similarly,   for  the  second  approximation  to  (E. 
terms  linear  in  C  and  O  only  are  required.     Higher  order  requires 
higher  powers. 

2.  5    Average  Temperature  Rise  Through  Reactor 

The  total  enthalpy  rate  of  increase,  AG,   in  the  exhausting  gases 
is: 

AQ  -  M  <2PCt;  ~r0) 

.(rv-Tj  z-rra*-  (cv  rv\   *|  A  >[  (2.44) 

V 

where  m  is  defined  as  the  total  mass  flow  rate  found  by  integrating 
the  individual  coolant  mass  flow  rates  over  the  cross  -  section  of  the 
reactor,   the  integrated  average  value  of  C     is  used,   and   Ou  is  the 
radius  of  the  reactor. 

However,    AQ  may  be  found  directly  ae: 


"   CP^ 


-18- 

o 

Therefore,   by  equating  Equations  2.  44  and  2.  45  we  may  solve  for  the 
average  temperature  rise  through  the  reactor: 


T-t.  ■ 


P  p.\/R?    i*^)  v^yy^ 


(2.46) 


The  exit  temperature  at  the  centerline  of  the  reactor  is  found 
by  putting    >1    equal  to  zero  in  Equation  2.  34: 

v+t.   -   <r„^fL     '       3. 

<x         \/fc(oj  (2.47) 


=       A(°^    T-    K°  (2.47*) 


Hence,  dividing  Equation  2.  47  by  2.  46,  the  ratio  between  the  center- 
line  exit  temperature  of  the  reactor  and  the  average  exit  temperature 
is: 

(2.48) 


T>)   -T.  ,Co)  f>j/fj^ 


T    -  T  <* 


/e(o7     ['{jf*})  i^y 


It  is  not  convenient  to  express  the  relation  between  the  maximum 
wall  temperature  and  the  average  temperature  increase  across  the 
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reactor  in  general.     However,   If    f(|)    and    g(  r]  )    are  known  functions, 
this  relation  can  be  found  quite  readily  as  will  be  shown  in  the  following 
examples. 


-20- 


m.     EXAMPLES 

3.  1    Reactor  Calculation  Using  Specific  Power  Distribution  Functions 
Thus  fax  in  this  thesis  the  general  equations  for  determining  the 
overall  temperature  rise  and  pressure  drop  through  the  reactor  have 
been  obtained  while  leaving  the  specific  power  distribution  as  general 
functions.     An  example  will  now  be  given  showing  how  the  results  ob- 
tained in  the  general  derivation  might  be  applied  to  a  given  reactor. 
The  general  solutions  for  the  neutron  density  distribution  in  a 
cylindrical  nuclear  reactor  are  given  in  terms  of  trigonometric  and 
Bessel  functions.     However,  for  convenience,  a  parabolic  neutron 
density  was  assumed  since  the  integrated  area  under  a  parabola  dif- 
fers but  slightly  from  the  more  generally  accepted  solutions  of  Bessel 
functions  or  trigonometric  functions.     For  this  example,   assume  a 
right- circular  cylindrical  reactor  having  parabolic  normalized  power 
distribution  characteristics  as: 


fC5)  =   '  -4e(5-i^  {yi) 

cj(Y|)     =      I    -4>V~  (3.2) 


Further,   assume  the  porosity  ie  uniform  across  the  reactor.     The 
power  distribution  parameters  to  be  perturbated  will  be   6    and    <$>, 
which  determine  the  relative  "flatness"  of  the  axial  and  radial  power 
input  cf  the  reactor.     These  parameters  are  primarily  a  function  of 
the  amount  of  reflector  material  used  in  the  reactor  design  if  homo- 
geneous fuel  loading  is  assumed. 
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With  these  functions,     B      is  found  from  Equation  2.  30  to  be: 


bQ-  |*(f)a§   =  i  -  % 


(3.3) 


Using  Equation  2.  48  and  working  out  the  integrals:,  indicated  to 
the  second  approximation  in   /*T*  and  '*ith: 

€(o)  =  /Ap/Po  JJJL  -^_  — A_  3.  A 


<»-*0[fr  +  ^•(■•wj] 


we  find,  after  simplification*   that: 


T,-T0    x 


»,     _To 


^+ii^-('*K)(»-rr^:) 


.{3.  5) 


To  find  the  maximum  wall  temperature,    which  will  occur  on  the 

centerline  of  the  reactor,   use  Equation  2. 14a  to  determine  the  value 

of    |    to  put  in  Equation  2.  14.     From  Equation  2.  14a,    T  occurs 

m 
when 


i 


and  hence,   from  Equation  2.  14: 


4© 


T~*T  '« 


XU>) 


U^      ?A'^     4fl       ,1/         ^J 


2       "'*» 


(3.6) 


and  using  Equation  2.  47a: 

r 


~i  *•*-  - 


'*\ 


I 


T.^-T, 


% 


(ffi-+i)f\//i-+4«      pJ  *"*<* 


-k-i 


(3.7) 
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Hence,   given  a  maximum  wall  temperature  limitation  in  the 
reactor,  the  relation  (as  a  function  of  0,  $,  f  ,  and  p.)  between  the 
average  temperature  increase  of  the  gas  flowing  through  the  reactor 
may  be  found  by  multiplying  Equations  3.  7  and  3.  5: 

t~„,-t„  ,         f  ft  >hr'^fa^ 


T, -T.  0-»/,)0-*0         A!.».i2£&(,+jJt)fr--^ 

5.  .V  ««•     4a        (3i^3(2 


S  (S«-     4©        (3/        3(3 


(3.8) 


It  is  simpler  to  leave  Equation  3.  8  in  this  form;  however*  for  given 

values  of  the  parameters  of  the  e  [tuition,  it  reduces  to  a  simple 

T  — T  T         —  T 

quadratic  to  determine  ~ — —  for  a  given 


If  all  wetted  surfaces  of  the  channel  are  heat  transfer  r  •  r faces, 
and  from  the  Reynolds  analogy,  it  is  reasonable  to  assume  that  ps3. 
for  fully  developed  turbulent  flow.     Furthermore,   it  is  assumed  that 
A  =  ©  for  comparison  purposes,  though  neither  of  these  assumptions 
are  requirements  of  the  analysis.     The  results  obtained  from  using 
Equation  3.  8  are  shown  in  Table  I.     This  table  show  s    the  maximum 
wall  temperature  required  to  obtain  a  temperature  increase  of 
1000°R  and  1500°R  for  an  initial  gas  temperature  of  1000°R  for  var- 
ious values  of  p  and  $.     The  table  was  prepared  on  the  assumption 
that  <t>  -  9  auid  P  =  £r  •     The  required  temperature  increases  were 
selected  quite  arbitrarily,   but  they  furnish  an  indication  of  the  trend 
in  maximum  wall  temperature  versus  average  temperature  increase 
through  the  reactor.     From  Table  I,   it  can  be  seen  that  the  maximum 
wall  temperature  required  for  a  given  temperature  increase  of  the 
working  gas  is  markedly  increased  as  the  power  distribution  varies 
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away  from  a  linear  relation  (  <*>  equals  zero)  or  when  the  coefficient 
of  heat,    f*  ,   gets  very  low. 

3.  1.  1    No -orifice  Condition 

The  static  pressure  loss  may  be  found  by  using  Equations  2.  44 
and  2.  45: 

YKcp(f,-T0)    =     AQ 

(3.9) 


or 

rra. 


But,  from  the  simplification  procedure  (not  shown)  used  in  deter- 
mining Equation  3.  5,  it  was  found  that 


K  - 


oC T.-T.      /         AP/pe 


ri-% 


Therefore,   Equation  3.  9  can  be  expressed  as: 


i-    ,-. — 


P-         \    f  T.        t         '^JUtt^J        pa-  (3.10) 

However,  -,       is  just  the  mass  flow  rate,  p  u    ,  per  channel, 

ex.  fro.1-  J  *o  o     r 

so  Equation  3.  iO  can  be  reduced  further  to: 


*f-[h  tJ^(w/v.jje 


(3.  10a) 


Figure  1  shows  the  variation  in  static  pressure  loss  with 
\  G.         (which  is  proportional  to  the  mass  flow)  for  various  values  of 
p-   for  an  assumed  temperature  rise  ratio  of  unity  across  the  reactor. 
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3.  1.  2   Artificial  Flow  Restriction  (Orifices) 

For  installed  orifices  scaled  so  that  the  mass  flow  in  the  outer 
channels  is  restricted  to  the  correct  value  to  give  a  uniform  tempera- 
ture rise  across  the  reactor,  the  average  exit  temperature  will  be 
identical  with  the  exit  temperature  from  the  center  channel  of  the  un- 
restricted reactor.     Therefore,  from  Equation  3.  7 

T~t7      fl-e/5)|(V  *H V    *6     /»'       3^        a       u\       (         ' 

If  the  orifices  restrict  the  flow  correctly,   the  pressure  decre- 
ment they  provide  will  make  all  channel  pressure  losses  identical 
with  the  pressure  loss  in  the  centerline  channel,  which  Is  unrestrict- 
ed.    Therefore,   the  first  approximation  to  the  pressure  loss  ratio 
can  be  determined  from  Equation  3.  4  b»r  assigning  a  value  to  VeTbJ 
under  given  reactor  conditions. 

The  other  change  that  will  result  from  artificially  restricting 
the  channel  entrances  will  be  the  total  mass  flow.     This  total  mass 
flow  can  be  determined  directly  using  Equation  3.  11.     Since  the 
average  temperature  of  the  gases  is  the  exit  temperature,  a  constant 
across  tiie  reactor,  by  using  Equation  2.  34  we  are  able  to  determine 
m    if  it  is  recalled  that  it  is  a  factor  in  the  parameter,  YiT    .     Hence, 
integrating,   the  total  mass  flow  is: 

(3.  12) 

is  the  exit  temperature  of  the  centerline  chan- 

\ 
nel,   Equation  3.  12  can  be  rewritten  by  using  Equation  2.  47  and  the 


nr\  -- 

CpT0 

i 

However, 

since  T, 

is  the  e 
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fact  that     m  =  -£  \f^  '■ 


m   -  **<*.''■■&:  f<=(o)(i-  fyz.) 


(3.  13) 


■*•  ^    Rocket  Example 

Using  the  assumed  reactor  power  distribution  characteristics 
of  Equations  3.  1  and  3.  2,  a  rocket  engine  using  molecular  hydrogen 

the  working  fluid  was  analyzed.     No  dissociation  effects  were  con- 
sidered.    The  fuel  was  assumed  to  be  liquified  in  the  propellant  tanks; 
however,  due  to  cooling  requirements  of  the  nozzle  and  reactor  reflec- 
tor,  the  hydrogen  was  assumed  to  have  reached  a  gaseous  state  before 
entering  the  pumps  of  the  "combustion  chamber". 

Account  was  taken  of  the  power  required  to  drive  the  pumps  by 
using  a  portion  of  the  gas  from  the  reactor  exit  to  drive  a  turbine 
which  furnished  power  for  the  pumps.     The  temperature  of  the  gases 
as  they  were  bled  from  the  reactor  exit  was  reduced  to  acceptable 
turbine  operating  temperatures  by  mixing  it  with  the  required  amount 
of  coolant  hydrogen  to  lower  the  turbine  gas  inlet  temperature  to  an 


assumed  value  of  2000  R. 


(i  urb/ne 


Schematic   Draw/ng  of  pocket 
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For  this  example,   define  the  following  parameters: 
f"Hp      mass  flow  rate  of  pump 
frw      mass  flow  rate  of  nozzle 

<?        bleed  fraction  to  drive  auxiliary  turbine,  therefore, 

rhp  =  rr\r  ( •  +  <n 

TTj^     reactor  exit-entry  total  pressure  ratio  *  piT  /P| 

TTm    assumed  exit -entry  total  pressure  ratio  due  to  mixing 
loss  from  exit  temperature  gradient  *  Pa-r/p 

^r\      nozzle  exit-entry  total  pressure  ratio  (assumed  1.  0) 


The  thrust,    F,    of  a  rocket  is  given  by  the  equation: 

F   =    rv^  u4      +    A4  ( f4  -pa.J 


(3.  14) 


where    p     is  the  ambient  pressure  and  subscript  4  refers  to  exit 

r  a 

station.     The  specific  impulse,    I  ,   is  found  by  the  equation: 

r  -      F  -  -  _^l     -t    **/fii*.  i       b 

1    -       mtf  $  j-f/At    ?+~?a-]  (3.15) 

Using  the  defined  parameters,   the  nozzle  mass  flow  rate,     m,  , 
can  be  found  using  the  continuity  equation  and  the  perfect  gas  relations 
for  a  choked  throat  as: 

_i-    =     (»,T  TTRTTm  /  r. 


t  V     '2T 

If  the  expansion  ratio,   P*t  fn    (  is  specified,  the  nozzle -exit 
ratio  can  be  determined  from  the  following  equation  if  isentropic 
flow  is  assumed  in  the  nozzle: 
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A4 


(3.  17) 


where  /p4   ~    "",  ^Xr/p4     •     Equation  3.  17  uniquely  defines  the 

area  ratio  and  subsequently  the  exit  Mach  number  from  the  relation: 

*±  =  -J-  i-^-fi  -  -^*<)y    ~°    (3.i8> 

At  W4  fti    I  2  J 


Kence,    MA    ia  now  known. 

Since   M4  has  now  been  determined,  the  exit  velocity,    u.  ,  can 


be  found  if  TA  Is  known. 


T<    -    X<r  ( 


lf"LM» 


+  J-^-N\ 


'4 


But'  \  ■  T1t  -  T.T  (  .  +  T-^J 

»-  T ,  T  J 

Therefore;  _-       _    -    /  (     ,  TX-T.V,        r.,   mM  (3.19) 


^•t.(«+^«;)(m  u-nC) 


where      -^-— — -    is  the  static  temperature  increase  found  from  Equa- 
tion 3.  8,    T.    is  as  turned  equal  to  T,      .   since  the  entering  Mach  num- 
ber will  be  very  small,   and  My  was  determined  from  Fig.   6,   Dailey 
and  Wood  (2),  modified  to  account  for  the  friction  in  the  coolant  pas- 
sages. 

Knowing  T,     ,     T.   is  now  known,   and  the  exit  velocity,  u.,   is: 


U4    =    a4M4-   V**'4    M*  (3.20) 
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The  total  pressure  ratio  across  the  reactor  is  found  in  a  similar 
manner  using  the  perfect  gas  relations:  s 


^i-^-(i-^)h^<y" 


(3.  21) 


where    —   is  the  static  pressure  ratio  across  the  reactor  as  deter- 

P. 

mined  from  Equation  3.  10a  by  assigning  various  values  of  <£.  .  After 
having  determined  ~  ,  the  static  pressure  of  the  exhaust  is  given  by 
the  relation:  _  JL_ 


-K^^,^*;)    "  (3'22) 

The  total  mass  flow  rate  that  must  be  pumped  by  the  pumps  is 
known  if   e.  is  assigned,   since  from  Equations  3.  10  and  3.  10a: 

fro.9-  p.     v  '  (3.  23) 

In  order  to  pump  the  required  total  amount  of  mass  flow,   the  pumps 
expend  power  which  was  assumed  to  be  supplied  by  the  turbine.     Equa- 
ting the  turbine  power  output  to  the  required  pumping  power: 

fir  ™P  <rP  ATturb  ngtur(o 


p^p  r^t  (3.24) 


Therefore,  assuming  that  the  mixture  of  bled  gas  and  coolant  gas 
enters  the  turbine  at  2000  R  and  exits  at  the  reactor  entrance  tem- 


perature, the  amount  of  bled  gas  from  the  reactor  exit  is  given  by 

'it  ™p     H*. 


<^  ™r  =    -"F^T     r,D  (3.25) 
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In  working  out  the  example,   a  "rubber"  engine  was  assumed, 
i.  c.  ,  the  rocket  was  assumed  to  be  correctly  designed  for  each  vari- 
ation in  the  parameters.     A.  graphitic  reactor  is  used  as  the  basis  for 
assuming  the  high  maximum  wall  temperature  (Ref.    3)  since  it  was 
visualized  that  a  rocket  engine  would  not  have  a  long  operating  life. 
No  account  was  taken  of  any  possible  mass  addition  to  the  flow  from 
the  reactor  material  itself. 

Assumed  conditions 

p.  3  50  atmospheres 

p     =  0.  1  atmospheres  (  53,  000  feet) 

T        =       5500°R 

m 

tr      =         0.  90  for  no- orifice  condition 
m 

1.  00  for  orifice  condition 

^     a  0.5 


Tl  * 

500°R 

FV*- 

50 

i    * 

1.35 

c    « 
p 

3.  8  Btu/lb  °R 

Using  these  conditions,   the  exit-throat  area  ratio  is  5.  53  and  the 
exit  Mach  number  is  3.  17. 

Figure  2  shows  the  specific  impulse  versus  ,6  obtained  from  the 
calculations  for  $  »  Q  m  0.  5.     A  linear  power  distribution  (<f>  ■  zero) 
is  also  shown  for  comparison.     No  attempt  is  made  to  distinguish 
between  the  values  of  mass  flow  used  at  each  ft,   since  their  effect  on 
the  specific  impulses  was  negligible,  although  their  values  do  affect 
the  pressure  loss  across  the  reactor.     From  these  results,  the  spe- 
cific impulse  appears  to  be  primarily  a  function  of  the  temperature 


rise  imparted  to  the  gas  in  passing  through  the  reactor.     Therefore, 
to  get  high  specific  impulses  in  a  rocket,   it  would  seem  preferable  to 
attempt  to  increase  the  temperature  of  the  coolant  gas  by  increasing 
the  heat  transfer  coefficient,   although  this  lowers  the  allowable  mass 
flow  rate  to  preclude  very  high  pressure  losses,   or  choked  flow,   in 
the  coolant  passages  due  to  the  higher  heat  input  to  the  gases.     For 
this  same  reason,    orificing  appears  to  be  advantageous  in  increasing 
the  specific  impulses. 

3.  3    Ramjet  Example 

In  the  ramjet  example,  the  conventional  burner  is  replaced  by  a 
nuclear  reactor  having  the  power  distribution  characteristics  of  Equa- 
tions 3.  1  and  3.  2.     This  analysis  was  based  on  standard  ICAO  condi- 
tions at  1000  feet  altitude. 

In  working  out  the  e;cample,   a  "rubber51  engine  was  assumed  and 
the  nozzle  was  assumed  to  be  designed  to  e.vpand  to  atmospheric  pres- 
sure unless  the  throat  static  pressure  was  lower  than  ambient.     In 
such  a  case,  no  expansion  was  considered,   and  the  nozzle  was  assumed 
to  be  terminated  at  the  throat. 
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For  this  analysis,   define  the  following  parameters: 

Tl\       ratio  of  free -stream  total  pressure  to  static  a  ffer/po 

TTj       ratio  of  total  pressure  across  diffuser  a    f**Vpo 

TTR       ratio  of  total  pressure  across  Reactor  a   f»T/pir 

^m      ratio  of  total  pressures  due  to  assumed  mixing  loss 
in  noa- restricted  case  a   pa'T/p        (Assumed  1  for 
restricted  case.  )  '  3r 

TTn     ratio  of  nozzle  total  pressures  (aasumed  1) 
The  thrust*    F ,  of  a  ramjet  is  given  by: 

Due  to  the-  fact  that  ramjet  engines  must  operate  for  long  peri- 
ods of  time,   and  in  an  oxidizing  atmosphere,   a  different  type  reactor 
material  must  be  used  in  a  ramjet  than  in  a  rocket.     Because  of  this, 
the  materials  suitable  for  ramjet  applications  can  withstand  a  much 
lower  maximum  wall  temperature  than  those  for  a  rocket  engine.     The 
heat  input  to  the  coolant  gases  will  be  much  less,  and  consequently  the 
acceleration  of  the  gases  in  the  coolant  passages  will  be  less.     There- 
fore,  in  this  example,   it  was  assumed  that  the  static  temperature 
rise  and  the  static  pressure  loss  across  the  reactor  were  approximately 
equal  to  their  total,   respectively. 

u<ra     can  be  found  by  finding  the  ratios  of  exit  Mach  number  to 
the  free  stream  Mach  number  and  exit  total  temperature  to  the  initial 
total  temperature. 


2. 

') 

(3.  27) 
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But:  v 

frr  -  M '  *  Hr  O  ^~  -  K  *d  ff«  ""..  ( '  +  *r  M°") "         <3-  28> 


and  from  this  we  obtain: 

be) =  (^  w)  ^F^  "  "5m?      (3-29) 


Also:  -T-  _ 

^  TOT    -  TOT  (3.30) 


AT     .    f,  -T 
where      ~^f~   -  — = —   of  the  general  reactor  analysis. 

Therefore,   the  exit  velocity  can  be  expressed  as: 


5  -/(' +  ^ V +  ¥<  L' ~  T^w^)^)1     <3-31' 


Since  it  was  assumed  that  the  static  pressure  los3  and  total 
pressure  loss  are  approximately  equal,    w~   is  found  by  using  Equa- 
tion 3.  10a  to  determine  Ap/p  for  various  values  of    £  . 

^R.    ~      l  ^piT  (3.32) 


Knowing      ~  ,   the  throat  pressure  can  be  found  from: 


h  '       =     M-    *"'   IT.  IQ   T.  T. 


/p«  =    777       "^   R  *•  (3-33> 


If  this  ratio  is  less  than  unity,   no  further  expansion  was  assumed; 

if  the  ratio  was  greater  titan  unity,   ideal  expansion  to  ambient  pressure 

was  assumed. 
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Figures  3  and  4  show  the  results  obtained  for  C_,  (based  on  the 

J? 

reactor  cross -sectional  area),   the  thrust  coefficient,   versus  Mach 
number  for  various  values  of  3  at  an  assumed  maximum  reactor  wall 
temperature  of  2960°R.     CF  is  defined  as  the  total  thrust  divided  by 
the  free  stream  dynamic  pressure  and  the  cross  -sectional  araa  of  the 
reactor.     Values  of  -rr «  used  in  the  example  were: 

t.S    2    S  o8s 

^.s  —  *.s-  o-  7£ 

rr     ,   the  assumed  mixing  loss,  was  arbitrarily  assigned  a  value  of 
0.  90  for  the  non- restricted  example  and  unity  for  the  orificed  exam- 
ple.    Results  shown  are  for  an  assumed  reactor  porosity  of  0.  5. 
The  curves  shown  are  for     €.  m  0.  03,   but  other  mass  flow 
rates  give  similar  curves.     The  one  shown  is  about  optimum,    since 
it  was  found  that  further  increase  in  mass  flow  rate  increased  the 
pressure  loss  across  the  reactor  to  a  point  where  the  overall  thrust 
is  reduced.    While  it  should  be  stressed  that  the  initial  conditions  and 
reactor  configuration  chosen  were  arbitrary  and  this  example  was 
worked  to  demonstrate  the  methods,    still,   the  trend  of  the  results 
should  be  realistic,   if  not  their  absolute  magnitudes.     From  these 
two  figures  at  the  arbitrarily  chosen  configuration,    optimum  per- 
formance is  obtained  at  p   of  approximately  3-5  at  a  Mach  number  of 
about  3.     Figure  4  shows  that  while  CF  at  the  optimum  Mach  number 
is  slightly  increased  by  orificing,   the  vehicle  performance  drops 
more  sharply  away  from  the  optimum  Mach  number,    so  it  is  doubtful 
that  there  is  much  advantage  to  orificing  for  a  ramjet. 
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TA3LE I 


Maximum  Wall  Temperature  Required  to  Obtain  Given  Average 

Temperature  Increase  of 
Coolant  Gas  (Initial  gas  temperature  assumed  to  be  1000   R) 


AT 

*  1000 

AT  *  1500      [ 

M*f 

Tx{o) 

T 
w 

m 

lyo) 

T 
w 
m 

20 

3430 

3440 

4870 

4830 

10 

3460 

3470 

4930 

4940 

<|>  3   1.  0 

5 

3520 

3625 

5040 

5250 

1 

3860 

7000 

5600 

10,630 

20 

2790 

2818 

3770 

3820 

10 

2801 

2868 

3810 

3910 

<>  =  .75 

5 

28/8 

2981 

3840 

4080 

1 

2950 

4760 

4030 

6850 

20 

2418 

2518 

3160 

3310 

10 

2420 

2533 

3170 

3340 

4  ■  .  50 

5 

2430 

2620 

3190 

3440 

1 

2480 

3720 

3260 

5170 

20 

2173 

2450 

2771 

3190 

10 

2175 

2452 

2772 

3190 

4>  a  .  25 

5 

2177 

2490 

2779 

3250 

1 

2193 

3190 

2801 

4315 
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TABLE    I 
(continued) 


AT  ■  1000 

AT  x  ] 

1500 

MPf 

T^o) 

T 
w 

m 

Tjfo) 

T 
w 

m 

20 

2000 

2000 

2500 

2500 

10 

2000 

2000 

2500 

2500 

4>  -  .  0 

5 

2000 

2000 

2500 

2500 

1 

2000 

2000 

2500 

2500 
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FIG.  I   -  STATIC     PRESSURE    LOSS    VS.  DIMENSIONLESS 

MASS    FLOW   RATE    FOR   UN  RESTRICTED  COOLANT 

PASSAGES 
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FIG.  2  -  SPECIFIC     IMPULSE     VS.   HEAT   COEFFICIENT    FOR 
NUCLEAR- HYDROGEN    ROCKET  EXAMPLE 
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FIG. 3  -  RAMJET     EXAMPLE  (UNRESTRICTED  COOLANT 
PASSAGES)  -  THRUST   COEFFICIENT   VS.  FREE 
STREAM    MACH    NUMBER 
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FIG. 4-  RAMJET  EXAMPLE  (ORIFICED  COOLANT 
PASSAGES  ) -THRUST  COEFFICIENT  VS. 
FREE    STREAM     MACH    NUMBER 


